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Wc present a formalism for self-calibrating tomography of arbitrary dimensional systems. Self- 
calibrating quantum state tomography was first introduced in the context of qubits, and allows the 
reconstruction of the density matrix of an unknown quantum state despite incomplete knowledge of 
the unitary operations used to change the measurement basis. We show how this can be generalized 
to qudits, i.e. d-level systems, and provide a specific example for a V-type three- level atomic system 
whose transition dipole moments are not known. We show that it is always possible to retrieve 
the unknown state and process parameters, except for a set of zero measure in the state-parameter 
space. 
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I. INTRODUCTION 

The characterization of quantum states is a key step in 
quantum information processing and quantum comput- 
ing [1]. Such characterization can be achieved by per- 
forming a variety of measurements on many identically- 
prepared copies of the state and using these results to 
reconstruct the quantum state. This process is known as 
quantum state tomography (QST) [2-4]. 

The various measurements are typically achieved by 
applying a unitary transformation to the quantum state 
before measurement in a convenient basis. This re- 
quires a well-characterized unitary operation, which is 
not always available. With such conditions in mind 
self- calibrating tomography (SCT) was introduced by 
Branczyk et al. [5], which allows the reconstruction of a 
quantum state despite incomplete knowledge of the uni- 
tary operations used to change the measurement basis. 
SCT was used to successfully reconstruct the state of a 
polarization-encoded photonic qubit using a wave plate 
with unknown retardance. In this work we formalize the 
ideas presented in [5] and show how to generalize SCT to 
the case of multi-level systems. 

The concept of performing tomography with unknown 
parameters in both the state, as well as the process, has 
also recently been investigated by Mogilevtsev et al. [6] 
and Merkel et al. [7]. 

II. TWO-LEVEL SYSTEM 

In this section, we formalize the ideas of self-calibrating 

tomography (SCT) in the context of qubits. Wc explore 
several scenarios including the one originally discussed 
in [5] and we prove that such schemes will be able to 
reconstruct any possible state of a qubit system except 
for a set of zero measure. In this exception, a different 
scheme can be used. 



A requirement of quantum state tomography is the 
ability to perform measurements in a variety of differ- 
ent bases. In many physical systems, there is a "pre- 
ferred" basis in which a measurement can be made, e.g. 
the horizontal/vertical polarization basis in polarization- 
encoded single photons; or the ground/excited state basis 
for atomic systems. Other bases are typically accessed by 
performing unitary operations on the state before mea- 
surement, effectively changing the measurement basis. 

Self-calibrating tomography considers the case where 
this unitary operation is not completely characterized, 
but controllable to some extent. It is "self-calibrating" 
in the sense that certain parameters of the unitary oper- 
ation do not need to be calibrated in advance. 

To demonstrate the concept, we begin by defining a 
general qubit density operator and a general generator of 
rotations that is related to the Hamiltonian that governs 
the unitary evolution of the state. We then show how 
the measurement statistics depend on the parameters of 
the initial state and the unitary. 

The single qubit state p and generator of rotations 
under which the state evolves G = j^J dtH{t), can be 
written in the eigenbasis of the Pauli operator ct^, i.e. 
{|0) ,|1)}, as follows: 

The components of the Bloch vectors associated with 
p and G are given by: 

Vp = {2poi cos(7), 2poi sin(7), poo - Pii} ; (2a) 

= {he CQs{(j)),hc sin{(j)),hz} . (2b) 

For p to provide a valid description of a quantum state, 
it must have positive eigenvalues and unit trace. Given 
the parametrization in Equation (la), these constraints 
are equivalent to 
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= V 4Poi + (POO - Pll) ^ < tr(p) = Poo + Pll = 1 . 

(3) 
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FIG. 1. The Bloch vectors vp and Vq, given in Equation (2), 
shown on: (a) the Bloch sphere; and (b) a cross-section of the 
Bloch sphere. The unitary operation Uv, shown in Equation 
(5), can be thought of as either: (c) rotating Vp about Vq\ 
or (d) rotating the reference frame about Vq in the opposite 
direction. 

The action of the Hamiltonian on the state will be 
given by 

p^p^= U^PqUI , (4) 

with 

Z^.(/i.,/ie,<^)=e-^<=(''-''-^). (5) 

We use the label v to refer to a specific set of parameters 
{hz, he, (/)} witin the generator G. 

By expanding the exponential it is easy to see that the 
Bloch vector of the state p is rotated by an angle 

n = \vq\ = ^W+Ki (6) 

around the axis specified by the Bloch vector of G. This is 
then an active transformation of the vector. Equivalently 
the action of the Hamiltonian can be defined in terms of 
the basis elements {|0) , |1)}: 

\Q)~^\0,)^Ul\0); (7a) 
\l)~,\U)^Ul\l). (7b) 

In this interpretation, the Bloch vector of the state re- 
mains fixed but the Bloch sphere itself rotates by an angle 
—Q, around the axis v^. This is the passive transforma- 
tion of the reference frame, which is equivalent to the 



transformation of the basis vectors and is illustrated in 
FIG. (1). Either way, transformation of p according to 
lAy, followed by measurement in the {|0) , |1)} basis con- 
stitutes measurements given by the projectors 

nl^ = |o,) (0,1 ; ni = |i,) (1,1 . (8) 

We draw attention to the fact that the superscript j in 
is a label denoting the original projector and should 
not be mistaken for an exponent. 

The statistics associated with measuring 11° are given 

by 

n°=tr(p'Il",) =^/,,p,,, (9) 
i<i 

where p' = Mp and A/" is a constant that depends on the 
duration of data collection, detector efficiency, loss etc. 
In what follows we reabsorb this constant into p' and thus 
recast the condition tr(p) = 1 into tr(/5') = N . Once the 
unnomarlized density matrix p' is determined it suffices 
to divide by its trace to recover p. The coefficients in 
Equation (9) are given by 

/oo = + s'hl ■ (10a) 
hi = s^hl ; (10b) 
/oi = 2 (csin(/3) + scos(/3)/i,) 7/^, (IGc) 

where hi = hi/il, c — cos{fl/2), s — sin(51/2), and /3 = 
(/) — 7. Similarly, the statistics associated with measuring 
II;', are given by 

< = tr(p'ni) =^5,,p,,, (11) 

where 

goo = s^hl ; (12a) 

ffn =c^ + s^hl ; (12b) 

goi = -2 (^csin(/3) + scos{/3)hz^ Vsoo ■ (12c) 

From the form of n[J and nj,, it can be seen that just 
one type of measurement, either 11° or 11;',, is sufficient to 
access all the elements of the density matrix in Equation 
(la). An appropriate set of different realizations of one 
of the projectors 11^ can be used to construct a set of 
equations invertible for the unknown parameters of the 
density matrix, thereby realizing quantum state tomog- 
raphy of the two-level system. 

Inspection of the coefficients /oi and goi shows that 
the system of equations can only be solved for /3. If the 
phase (j) is known, then the system is determined and a 
complete inversion can be performed. 

We now define a set of measurement statistics A C 
{n1_^ , . . . , , n},^ . . .} and a set of unknown param- 
eters r C {poo, Poi, Pii,hz,hc, P}. Determining whether 
the set of equations in (9) and/or (11), parameterized 
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by A, will be invertible for F is a standard problem of 
analysis [8]. This occurs when the Jacobian, i.e. the 
determinant of the Jacobian matrix, is non-zero, 



J = det 



OA 



^0. 



(13) 



In the remainder of this section, we discuss several sce- 
narios, beginning with a completely characterized unitary 
before extending the formalism to include unknown pa- 
rameters in the generator G. 



A. Completely characterized unitary 

We first consider the familiar scenario where one has 
complete knowledge and control over the unitary oper- 
ation. Given the ability to measure in one particular 
basis, one can make measurements in any arbitrary basis 

{|o.),IM}. 

As an illustrative example, we consider the four mea- 
surements introduced by White et al. [9] which corre- 
spond to the four measurement operators 



n° = |0) (0| ; 
ni= |1)(1| ; 

ILl=U^\l){l\Ul; fori. = i, i 

where the unitary operations are given by 

Ui=U{0,Tr/2,0); 
Z^ii=ZY(0,7r/2,7r/2), 



(14a) 
(14b) 
(14c) 

(15a) 
(15b) 



and i4{hz-i he, (j)) is defined in Equation (5). The Jacobian 
for this scenario is given by 



Ja = det 



a(n°,nSniSni\) 
5(Poo,Poi,Pii,7) 



Poi • 



(16) 



As long as p has non-zero coherence, we can solve for all 
unknown parameters F = {Pooj Poi) Pii) 7}- H Poi = 0, 
the phase 7 will be undefined, however, the system of 
equations will still be invertible as in this situation 7 will 
be absent from the expressions for n^. 



B. One unknown parameter 

We now consider a scenario where G contains one un- 
known parameter. To solve for this additional parameter, 
we require one additional measurement. We consider the 
scenario introduced by Brahczyk et al. [5] (with some mi- 
nor adjustments), which requires hz = 0, and he = rricXc 
where nic is a controllable parameter and Ac is unknown. 
We also require the ability to control 0. Such a physi- 
cal scenario can be realized when a polarization encoded 
qubit passes through ttIc wave-plates of unknown retar- 
dance Ac at an orientation specified by (j). 



Under these conditions, we construct five measurement 

operators: 



ni= |1)(1| ; 

Ul =U^\l) {l\Ul ; for 1/ = i, ii, iii, iv 
where the unitary operations are given by 

Ai=Z^(0,2Ac,0); 

Aii=zi(0,Ac,V2); 

=W(0,2Ac,7r/2), 

and Uihz, he, (f) is defined in Equation (5). 
The Jacobian for this scenario is given by 



(17a) 
(17b) 

(18a) 
(18b) 
(18c) 
(18d) 



det 



/ d{n^,nl,nl^,n\^^,nl^) 
\5(poo,Poi,Pii, Ac,7) 

64poi sm'^ \ — \ cos* 
X (sin(7) — cos(7)) . 



(19) 



(20) 



Jb will be equal to zero in the following circumstances: 

Poi=0; (21a) 

Ac = 0; (21b) 

Ac=7r/2; (21c) 

7 = 7r/4, (21d) 

and therefore the transformation will not be invertible. 
Let's examine each circumstance individually. The case 
where poi = was discussed in the previous section. 
When Ac = 0, we never change the measurement basis 
and thus trivially nothing can be learnt about the state 
except for the value of poo- If -'^c = ti": then lA^^ = Ui^, lead- 
ing to redundancy in the measurement statistics thus pro- 
viding insufficient information. To understand the last 
condition, one must look at the explicit forms of nl and 
n?;, which differ from n?; and n;\, by sin(7) cos(7) re- 
spectively. When 7 = 7t/4, sin(7) = cos(7), once again 
leading to redundancy in the measurement statistics. 

In all other circumstances, the system will be solvable 
for the unknown parameters in the density matrix and 
the unitary. 



C. Two unknown parameters 

Wc now introduce an additional unknown parameter 
in the generator, such that h^ = nizXz where ruz is con- 
trollable and A^ is unknown. We can build upon the 
results of the previous section, where we had implicitly 
set niz = 0, and make an additional measurement given 
by the measurement operator 



fll=U^\l){l\Ul 



(22) 
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FIG. 2. Representation of two qudit systems with d = 3. 
In (a) each level can be coupled to each other level by the 
Hamiltonian interactions represented by the blue arrows. In 
(b) the two levels |1) and |2) cannot be coupled directly but 
can still interact to second order through the level |0). 

where the unitary operation is given by 

U^^U{K,Q,Q). (23) 
The Jacobian for this scenario is given by 

Jc = ^Jb = 6A,VooJb, (24) 

where Jb is defined in Equation (20). The Jacobian 
adopts this simple form because Jc is upper block trian- 
gular where the first block is identical to Jb and the other 
block is simply given by dn\/dXz- As long as Jb ^ 0, and 
Az ^ in the sixth measurement, the system will be solv- 
able for the unknown parameters in the density matrix 
and the unitary. 



can be broken down into performing a succession of two- 
qubit SCT on all possible pairs of levels in the system, 
denoted by ovals in FIG. 2(a). 

On the other hand it might happen that there are 
eigenstates that are not directly coupled, but are coupled 
to another common state. This is the case illustrated in 
FIG. 2(b). In this scenario, pairwise qubit tomography 
will not be sufficient. Nevertheless it is still possible to 
gain complete information about the state and, in partic- 
ular, the value of the coherences between the levels that 
are not directly coupled, by employing ancillary levels 
such as the |0) state in FIG. 2(b). This will be demon- 
strated in the next section, for a three-level atomic sys- 
tem in the V configuration. 

We note that the general idea of how to perform quan- 
tum state topography of d-level systems have been stud- 
ied by Thew et al. [19]. 

IV. EXAMPLE: V-TYPE THREE-LEVEL 
SYSTEM 

In this section we provide an example of how to extend 
SCT to the case of a three-level system in which two of 
the levels are not directly coupled by the Hamiltonian 
that is used to perform changes in the measurement ba- 
sis; FIG 2 (b). We consider a concrete physical system, 
nevertheless the conclusions reached in the following sub- 
sections will be independent of the implementations that 
we outline here. 



III. GENERALIZATION TO 
MULTI-DIMENSIONAL SYSTEMS 

The generalization to multi-dimensional systems fol- 
lows readily. One applies the same procedure: 1) write 
down the unitary being applied to the state, making note 
of parameters that are unknown as well as parameters 
that are controllable; 2) construct a set of measurement 
operators analogous to those in Equations (17) using dif- 
ferent realizations of the unitary, achieved with differ- 
ent settings of the controllable parameters; 3) check that 
the measurement operators are suitable by calculating 
the Jacobian in Equation (13); 4) collect measurement 
statistics from measurements given by the measurement 
operators in step 2; 5) solve for the unknown parameters 
using one of a variety of methods, including simple linear 
inversion, least-squares estimation or the popular maxi- 
mum likelihood estimation method [10], or alternatively, 
one of a growing number of new techniques [11-18]. 

The above procedure may be simplified if the multi- 
dimensional system supports interaction Hamiltonians 
that generate all combinations of pair-wise coupling be- 
tween the eigenstates, as illustrated for a three-level 
atomic system in FIG. 2(a). In this situation, each 
measurement operator can be constructed from a two- 
dimensional unitary operation, and the SCT procedure 



A. Description of the system 

We consider an atomic or molecular three-level sys- 
tem with a non-degenerate V-type energy ladder with 
excited states ]1) and ]2) and ground state ]0). The mat- 
ter Hamiltonian is given by: 

Ho = huji\l){l\+hu:2\2){2\ , (25) 

where the energy of the ground state ]0) has been set to 
zero. The system will couple to the electromagnetic field 
via the transition dipole moment operator d, where 

(26a) 
(26b) 
(26c) 

The light-matter interaction will be given in the dipole 
approximation by 



dii 


= ^22 


= di2 =^21=0 


doi 


= d*io 


^0; 


do2 


= ^20 


^0. 



(27) 



where E{t) = E^{t) + E~ [t) is the electric field. In par- 
ticular, it is assumed that the field is in a superposition of 
two continuous-wave lasers that are each resonant with 
the ]0) — > Jl) and JO) — >■ [2) transitions: 



E+[t) 



Eie''^'e'''^''ui 



(28) 
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where Ei and 9i are the controUable amphtudcs and 
phases of the fields respectively, Ui are their polarizations, 
and uji their frequencies. Without loss of generality, we 
take Ei to be real. In the "bare" basis of the V-type 
system, the light-matter interaction Hamiltonian is 

HUt) = d- E{t) (29) 
= (dio |1) (0| + ^20 |2) (0| + h.c.) • E{t). (30) 

Making the rotating wave approximation and going into 
the interaction picture with respect to the Hamiltonian 
Hq, the interaction Hamiltonian is given by 

H,{t) =(E^e'^^dw ■ ui + E2e"^'e'^'''d,o ■ U2) |1) (0| 

+ (^Eie*^'e'^''''d2o ■ ui + i^2e'^^(?20 • ^2) |2) (0| 



B. Measurement statistics 

We now calculate the measurement statistics associ- 
ated with projection onto the two excited states of the 
system. 

A simple way of parametrizing the density matrix of a 
3 level system is as follows: 



Poo 



Poie-'^"' P02e-''-'«' 



P 



Poie 



if 01 



Pii 



Pi2e 



-2712 



(38) 



P22 



+ h.c. 



(31) 



where A jj = Wi — uij . We now calculate the generator of 
rotations by assuming that the field is turned on at time 
—t/2 and turned off at time t/2: 



Note that unlike the qubit case, the constraints that 
the set of numbers pij will have to satisfy so that p is 
a non-negative linear operator are by no means trivial 
[21, 22]. 

To perform state tomography on the unknown state p, 
we evolve it according to Ui, = exp(— iG), where we as- 
sume that the interaction time is short compared with the 
rate of spontaneous emission. After the evolution is com- 
plete, the detection of a spontanously-emitted photon of 
frequency oji or UI2 constitutes a projective measurement 
of p given by the projectors 



. 1 

G=- dt'H,{t') 

" J -t/2 



11.) (i„ 



|2.) (2.1 



(39) 



(32) 



= Eite'^^dw ■ Ml |1) (0| + E2te'^^d2o ■ U2 |2) (0| 
+ h.c. (33) 

where we also assumed that t ^ such that 



where = Ui, The statistics associated with mea- 
suring II;', are given by 

2 

nl=tv(^p'Ul)=J2fijPij, (40) 



*/2 iAi2t 
■t/2 ^ 



/- 

generator as 



2fsinc(Ai2t/2) « 0. We can write the ^^^re p' =J\fp and the coefficients are given by 



G = y e*^i 1 1) (0| + y e^-^^ |2) (0| + h.c. (34) 



where 



hj = 2Ejt\djo ■ Uj\; 
'Pj = + arg(d^o ■ Uj) . 



(35a) 
(35b) 



Equation (34) can also be written in the ordered basis 
{|0),|l),|2)}as 



where 



^ ^ e''^2/i2 




(36) 



The matrix above has both zero trace and zero de- 
terminant, and therefore the eigenvalues take the very 
simple form spec(G) = {- + /ii/2, 0, ^hl + hl/2]. 
Another elegant way of writing the above matrix is in 
terms of the Gell-Mann matrices [20] : 



joo = «is ; 


(41a) 


fii = {hjc + hiy ; 


(41b) 


f22 = hlJ4{c-lf; 


(41c) 


foi = 2 sin (/3oi) \/ /oo/ii ; 


(41d) 


/02 = 2 sin (^02) V /00/22 ; 


(41e) 


/12 = 2 cos {(3i2) \J /11/22 , 


(41f) 


hj = hj/Q.] 


(42a) 


c = cos(17/2) 


(42b) 


s = sin(0/2) 


(42c) 


0= ^Jhl + hl 


(42d) 


Poj = <Pj - 7o,j ; 


(42e) 


P12 = 7l,2 + (I>1- (l>2- 


(42f) 



G = ^ (cos{(f)i)hiXi -|-sin((j!)i)/iiA2 

-I- C0s(^2)/l2A4 -I- Sm{(l)2)h2X5^ . 



Similarly, the statistics associated with measuring are 
given by 



(37) 



tr 



i<j 



(43) 
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where 



500 


= hls^ ; 


(44a) 




1 z V / ' 




922 


= {h\ + hl^'' ; 


(44c) 


901 


= 2 sin (/3oi)\/5oo5ii; 


(44d) 


902 


= 2sin (^02) Vfl-ooto ; 


(44e) 


912 


= 2 cos (/3i2) s/ 911922 ■ 


(44f) 



Prom the form of nl and n^, it can be seen that just 
one type of measurement, either 11;', or 11^, is sufficient to 
access all the elements of the density matrix in Equation 
(38). Different realizations of one of the projectors 11^ 
can be used to construct a set of equations invertible for 
the unknown parameters of the density matrix, thereby 
realizing quantum state tomography of the three-level 
system. 

We note that, as in the qubit case, the equations are 
only invertible for (3ij , but not for the individual quanti- 
ties jij and (j)k ■ If the phases 4>i and (/>2 are known then 
the system is determined and a complete inversion can 
be performed. 

Physically, both fields which couple |1) and |2) to the 
ground state must be present simultaneously in order to 
access all elements of the density matrix. This is to be ex- 
pected since the only way to have a coherent interaction 
between the two excited states in a V-type system is by 
going through the ground state — possible only when both 
fields are present. This can also be seen from the equa- 
tions above, since /12 and gi2 are dependent on /11/22 
and 91x922 respectively. 

On a more abstract level, it is not possible to gener- 
ate the whole set of Gell-Mann matrices from the sets 
{Ai,A2} or {A4,A5}, however the entire set can be ob- 
tained by multiplying pairs from the set {Ai, A2, A4, A5}, 
which is precisely the set contained in the expression for 
G in Equation (37). 



C. Unknown transition dipole moments 

In this section, we demonstrate how self-calibrating 
tomography can be performed on a V-type three-level 
system with unknown transition dipole moments djQ. 
Without loss of generality, wc assume that the transition 
dipole moments are real and positive. We can therefore 
parameterize the generator in Equation (34) as follows: 

hj — rrijXj ; (45a) 
<t>j = Oj , (45b) 

where nij = 2Ejt is a controllable parameter propor- 
tional to the square of the intensity of the fields, Xj = 
\djo • Uj\ is unknown and 6j is the controllable phase of 
each field. 



To demonstrate the process of constructing the appro- 
priate measurement operators, we break the problem up 
into three parts. The first part consists of treating one 
ground-excited state pair as a two-level system; this gives 
5 measurement operators analogous to the qubit case dis- 
cussed above. The second part consists of doing the same 
for the other ground-excited state pair; this gives four 
additional measurement operators since the ground state 
population parameter is already obtained from the first 
part. The third part consists of cleverly picking two more 
operators to solve for the coherences between the two ex- 
cited states. 

The first set of measurement operators are given by 



ni=|l)(l|; (46a) 

III =U^\1} {l\Ul; for 1/ = i, n, iii, iv , (46b) 

where the unitary operations are given by 

Wi =W(Ai, 0,0,0) ; (47a) 

=W(Ai,0,0,7r/2); (47b) 

Aii =ri'(2Ai, 0,0,0); (47c) 

Ziiv =Zi(2Ai,0,0,7r/2), (47d) 

where 

U{h,,h2, 01, 02) = c-'<=(''i-'^^>i>^) , (48) 



and G(/ii, /i25 01) ^2) is defined in Equation (34). The 
second set of measurements is given by 

til = |2) (2| Ul ; for 1/ = v, vi, vii, viii , (49) 

where 

Z^v =^(0,A2,0,0) ; (50a) 

Wvi =Z^(0,A2,0,V2); (50b) 

Wvii =W(0,2A2,0,0); (50c) 

Z^viii =W(0,2A2,0,7r/2). (50d) 

The third set will need to provide information about 
the coherences between the two excited states. For this, 
we require both fields to be present simultaneously. 

ill=U,\l){l\Ul; fori/ = ix, X, (51) 

where 

Zii, =W(Ai,A2,0,0); (52a) 

= W(Ai,A2,0,0 + 7r/2) . (52b) 

To verify that the measurement statistics = 
tr(n^/5), given by the above operators, generate a set 
of equations invertible for the unknown parameters, we 
calculate the Jacobian defined in Equation (13). We 
note that the three sets of inversions are independent 
of each other. Therefore, as long as one can retrieve: 
{PoO) Pii, Poi) Ai, 71} from the first set of measurements 
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Vv ^c.T' >1' ^n'V 




FIG. 3. Structure of the Jacobian matrix for the V-type three- 
level system. White squares indicate a null element of the ma- 
trix. Black lines denote the three blocks associated with Ji , J2 
and J3. Gray non-zero elements outside the black boundaries 
signify that each block can only be solved if the blocks before 
it have already been solved. 



and Ai — Ac. We note that in all cases, the Jacobians are 
proportional to the coherences simply because the phases 
7i J are undefined if the coherences pij — 0. 

We also note that the measurement operators pre- 
sented here are by no means unique. Other unitaries, 
such as those that always contain both fields, may be 
preferable from an experimental perspective. The ad- 
vantage of the above choice of operators is two-fold. 
Firstly, it demonstrates that the task of constructing 
measurement operators for larger-dimensional systems 
can largely be broken up into the less daunting task of 
performing SCT on two-level systems. Secondly, it allows 
for a more efficient computation of the Jacobian to test 
the choice of operators. 

This analysis can be extended to systems with more 
complex coupling structures. One can use perturbation 
theory to show that if two levels are coupled by d ancillary 
levels, an expansion to order d -I- 1 of the propagator will 
be able to extract the coherence between those levels. For 
the case of rf = 1 this reduces to the qutrit case discussed 
here. 



in Equations (46); {p22, P02, A2, 72} from the second set 
in Equations (49); and {^12,712} from the third set in 
Equations (51) the inversion is possible. 

Equivalently, one can compute the Jacobian matrix 
and notice that it is upper block triangular where the 
first 5x5 block corresponds to (46), the second 4x4 
block corresponds to (49), and the final 2x2 block cor- 
responds to (51). This is demonstrated schematically in 
FIG. 3. Explicitly the Jacobians for the blocks are: 

X (sin (701) - cos (701)) ; 
= 2sin^(A2)cos(A2)po2 (cos(7o2) - sin(7o2)) ; (53b) 

•[^gp^2 \2 „:„4 „ /'\2 , \2„„„/'i2\\2 



,^Aisin^(g)pi2 (Af + Aicos(f))^ 

IF ■ ^^'^> 

Given the block structure of the matrix, the Jacobian 
of the whole transformation is simply 



J — J1J2J3 ■ 



(54) 



The analysis of the singularities of the quantities above 
(i.e. when Ji = 0) is completely analogous to that for the 
qubit case. In fact, Ji = Jb in Equation (20) if 701 7 



V. SUMMARY AND CONCLUDING REMARKS 

We presented a general formalism for performing 
self-calibrating quantum state tomography of multi- 
dimensional systems. This technique is applicable in sit- 
uations where the unitary operations used to change the 
measurement basis are not completely characterized. 

We introduced a condition, based on the Jacobian ma- 
trix relating the measurement statistics and the unknown 
parameters in both the state and process, to ensure that 
the choice of measurement operators will lead to a set of 
invertible equations. 

We then showed how this technique can be used to per- 
form quantum state tomography of a V-type three-level 
system with unknown transition dipole moments. In do- 
ing so, we demonstrated that the task of selecting ap- 
propriate measurement operators for multi-dimensional 
systems can be broken down into the less daunting task 
of choosing measurement operators for qubit SCT. 
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